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The relation between the on-sheU quark mass and the mass defined in the modified minimal 
subtraction scheme is computed up to order aj,. Implications for the numerical values of the top 
and bottom quark masses are discussed. We show that the new three- loop correction significantly 
reduces the theoretical uncertainty in the determination of the quark masses. 
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In Quantum Chromodynamics (QCD) practical calculations are very often performed in the modified minimal 
subtraction (MS) scheme [^|| leading to the definition of the so-called short-distance MS mass. The MS mass 
^ ' occupies a distinguished place among various mass definitions. First, it is a truly short distance mass not suffering 
On • from nonperturbative ambiguities. Second, the MS mass proves to be extremely convenient in multi-loop calculations 
' of mass-dependent inclusive physical observables dominated by short distances (for a review see ||^). On the other 
hand the experiments often provide masses which are tightly connected to the on-shell definition. Thus, conversion 
formulae are needed in order to make contact between theory and experiment. The two-loop relation between MS and 
^\ , the on-shell definition of the quark mass has been obtained in Until recently the accuracy of this equation was 
On ' enough for the practical applications. Meanwhile, however, new computations have become available which require 
, the relation between the MS and on-shell mass at 0(q!^) in order to perform a consistent analysis. The necessity of 
O i' an accurate determination of the quark masses, especially those of the top and bottom ones, is demonstrated by the 
' following two examples: 

(i) The main goal of the future B physics experiments is the determination of the Cabibbo-Kobayashi-Maskawa 
f~| , matrix elements which will give deeper insight into the origin of CP violation and possibly also provides hints to new 
physics. In particular the precise measurement of Vcb is very promising. It is determined from semileptonic B meson 
• ^ , decay rates. Thus it is desirable to know the bottom quark mass as accurately as possible as it enters already the 
' Born result to the fifth power. 

{a) One of the primary goals of a future electron-positron linear collider (NLC) or muon collider (PMC) will be 
the precise determination of the top quark properties, especially its mass, Mj. In hadron colliders like the Fermilab 
TEVATRON or the Large Hadron Collider (LHC) the top quarks are reconstructed from the invariant mass of the W 
bosons and the bottom quarks. On the contrary in lepton colliders it is possible to determine the top quark mass from 
the line shape of the production cross section cr(e+e~ — > it) close to the threshold. Simulation studies have shown 
that an experimental uncertainty of 100 — 200 MeV in the top mass determination can be achieved Thus also 
from the theoretical side the ambiguities have to be controlled with the same precision. In particular in this context 
much attention has been devoted to the relation of the pole mass. A/, to the MS mass, m. Although the pole mass 
demonstrates a bad infra-red behaviour it is often convenient to use it in intermediate steps. 
The connection between the MS and on-shell mass is given by 

m,{fi) = z,n{fi)M , (1) 
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where Zm is finite and has an expUcit dependence on the renormalization scale /i. In Q the infra-red finiteness and 
the gauge invariance of M was proven. In Q the perturbative expansion of Zm has been computed up to order a^. 
The main purpose of this letter is the computation of z„i up to order a^. Therefore three-loop corrections to the 
fermion propagator have to be considered. It is convenient to parameterize them in the following form 



i{M) = 1 + umCfz^ + all 



ry2 FF i n n FA . rp FL . rp FH 



'■M 



4')+0(al,), (2) 



with qm — a^r^' {M)/t:. ni is the number of light (massless) quarks. In the following wc will use z„i — Zm{M). 

The relation between the MS and on-shell mass is obtained from the requirement that the inverse fermion propagator 
has a zero at the position of the on-shell mass. Thus in principle it is necessary to evaluate three- loop on-shell integrals 
in order to obtain the 0{al) term in the MS-on-shell mass relation. This is avoided by computing expansions of the 
quark self energy for small and large external momentum. After building the proper combinations needed for the 
relation between the MS and on-shell mass a conformal mapping g^/M^ = 4a;/(l -I- w)^ is performed which maps 
the complex g^-plane into the interior of the unit circle in the w-plane. The relevant point lo — 1 (corresponding to 
~ M^) is reached with the help of Pade approximation |^]. For further details we refer the reader to Note that 
this method has already been extremely successful in other applications ||^. 

Let us at this point discuss the one- and two-loop results in order to get some feeling about the quality of our 
procedure. Using our method the following numbers can be extracted 

z^ = -1.0005(6), z^^ = -0.49(2), z^^ = -3.4(1), z^"^ = 1.566(4), z^,« = -0.1553(2), (3) 

where the error is obtained by doubling the spread of the different Pade approximants. The comparison with the 
exact result results {—1, —0.51056, —3.33026, 1.56205, —0.15535} shows very good agreement. 

The results obtained at one- and two-loop order encourage to apply the same procedure also at order a^. For the 
technical details we also refer to We just want to mention that in total we were able to evaluate 14 input terms 
for the Pade procedure. 

We could apply the described procedure to each colour factor separately For most practical applications in 
QCD the number of flavours may be set to 3, i.e. Cf = 4/3 and Ca — 3, and T = 1/2 can be adopted. Adding in a 
first step the results for the moments and performing the Pade approximations afterwards leads us to 

^^^^^ = 1 - 1.333aM + aij [-14.33 + 1.041n,] + al^ [-202(5) -I- 27.3(7)n, - 0.653n?] , (4) 
^ = 1 - 1.333aM + alt [-11.67 + 1.041n;] + al; [-170(5) + 24.8(7)n, - OMSnf] , (5) 

— ^ = 1 + 1.333a„. + al, [13.44 - 1.041n;] + af„ [l94(5) - 27.0(7)n; + 0.65371?] , (6) 
m[m) 

with a„i ~ ai"'^\m)/TT. For simplicity ^ — M and ^ — m has been chosen in (^) and (^, respectively. In the above 
equations the exact values of the nf term [^ is displayed. Our method leads to 0.656(8) which is in very good 
agreement. This is a further justification for our approach. 

The results for the values of n/ = 0, . . . , 5 are summarized in Tab. | where also the two-loop coefficients are displayed. 
Actually the coefficients of the terms linear in ni have been obtained by performing a fit to the three-loop results 
of Tab. |. The errors of about 2-3% for the three-loop results of (^)-(^) and Tab. | have again been obtained by 
doubling the spread of Pade approximants. On one side this is justified with the behaviour at 0{al) where the results 
of the first column in Tab. ^ are reproduced with an accuracy below 2%. On the other side the Pade approximants 
demonstrate more stability in the case where the moments are added and n; is fixed afterwards than in the case where 
the Pade procedure is applied to the individual colour structures separately. We want to stress that the errors in 
Eqs. (^)-(||) are somewhat overestimated once a value for rt; is chosen. Thus for practical applications Tab. || should 
be used. 

Let us next compare our results with various predictions which already exist in the literature obtained with the 
help of different optimization procedures. In Tab. ^ the results obtained for M/m(m) using the fastest apparent 
convergence (FAC) [|l^ and the principle of minimal sensitivity (PMS) [jl^ are compared with ours. For ni = 2 the 
discrepancy with our central value amounts to only 7%. It even reduces to 2% for n; = 5, i.e. in the case of the top 
quark. The predictions obtained in the large-/3o limit, where f3o is the first coefficient of the QCD (3 function, are also 
shown in Tab. Excellent agreement below 1% is found for ni = 3. It amounts to roughly 5% for ni — 4 and 14% 
for ni = 5. 
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In the remaining part of this letter we will discuss some important applications of the new term of O(a^) in the 
MS-on-shell relation. 

Threshold phenomena are conveniently expressed in terms of the pole mass. To be specific let us consider the 
production of top quarks in e'^e" collisions. The corresponding physical observables expressed in terms of Alt show 
in general a bad convergence behaviour. In the case of the total cross section the next-to-next-to-leading order 
corrections p^-p^ partly exceed the next-to-leading ones. Furthermore the peak position which is the most striking 
feature of the total cross section and from which finally the mass value can be extracted depends strongly on the 
number of terms one includes into the analysis. The reason for this is that the pole mass is sensitive to long-distance 
effects which results in intrinsic uncertainties of order Aqcd 0]- 

Several strategies have been proposed to circumvent this problem . They are based on the observation that 

the same kind of ambiguities also appear in the static quark potential, V{r). In the combination 2Mt + V{r), however, 
the infra-red sensitivity drops out. Thus a definition of a short-distance mass extracted from threshold quantities 
should be possible. The relation of the new mass parameter to the pole mass is used in order to re-parameterize the 
threshold phenomena. On the other hand a relation of the new quark mass to the MS mass must be established as it 
is commonly used for the parameterization of those quantities which are not related to the threshold. In order to do 
this consistently the three-loop relation between the MS and the on-shell mass is needed. 

In the concept of the so-called potential mass, mt,ps, has been introduced. Its connection to the pole mass 
is given by rrit^psifJ-f) — ^'^t ~ 5mt{^f) where Smt{fJ.f) is obtained from the static quark potential. In this way a 
subtracted potential, V{r,fj,f), is defined. The factorization scale ^/ has been introduced in order to extract the 
infra-red behaviour arising from the potential. In the combination Mt — SmtUjj) it cancels against the one of Mt 
leading to an significant reduction of the long-distance uncertainties in mt^ps fl^]- Thus it is promising to formulate 
the threshold problems in terms of rrit^psifJ-f) ^^'^ ^(^: M/) instead of Mt and V{r). In the final result the dependence 
on cancels. For the numerical analyses the value = 20 GeV has been adopted in [Q as its upper bound is 
roughly given by MtCpcisipi)- Smt{^f) is known up to order jlTj. Thus we are now in the position to establish a 
relation between the two short-distance masses mt.psipf) ^-nd mt{fJ.) with the result: 

mt,ps(20 GcV) = (165.0 + 6.7 + 1.2 + 0.28) GeV , (7) 

where the different terms represent the contributions of order to a^. For the numerical values rat{mt) = 165.0 GeV 
and ai^' (rrif (mf )) — 0.1085 have been used. Note that the error of the 0{al) coefficient in the MS-on-shell mass 
relation is negligible. The comparison of Eq. (|^) with the analogous expansion for Mi, 

Mt = (165.0 + 7.6+ 1.6 + 0.51) GeV, (8) 

shows that the potential mass can be more accurately related to the MS mass than Mt- The last term of the expansion 
in (0) is of the same order of magnitude as the error in the top quark mass determination at a NLG. Whereas in jlj] 
this term has been taken as uncertainty in the mass relation the error reduces significantly after the knowledge of 
the 0{al) term of the MS-on-shell relation. This can be deduced from the well-behaved expansion in Eq. (|^). The 
dominant error is now provided by the uncertainty in as- 

A similar strategy has been proposed in [|5|. There the so-called IS mass, M}^ , has been defined as half the 
perturbative mass of a fictions toponium \^Si ground state which would exist if the top quark were stable. The 
philosophy is very similar as in the case of the potential mass. From the experiment the quantity M^^ is extracted. 
In ||l^ it has been shown that this is possible with an uncertainty of approximately 200 MeV. In a next step M}^ has 
to be related to the MS mass mt{mt). As the extraction of M^^ is based on a next-to-next-to-leading order formalism 
the 0{a'l) relation computed in this work is necessary. 

In the practical calculation care has to be taken in connection to the expansion parameter which has to be used. 
Actually the so-called T-expansion has to be adopted. Details can be found in |l9y20[| . One finally arrives at the 
following relation between the MS and IS" mass 

TOt(TOt) = (175.00 - 7.60 - 0.97 - 0.14) GeV , (9) 

where Mt'^ — 175 GeV and ai^\Mz) = 0.118 has been adopted. Using the large-/3o results for the order term 
the last term reads —0.23 ||l^ which is off by more than 50% form the exact result. The conclusions which can 
be drawn from Eq. (^) are very similar to the ones stated above: the uncertainties due to unknown terms in the 
mass relations are negligible as compared to the error with which M^^ can be extracted from the experiment. The 
dominant uncertainty comes from the error in which amounts for ±0.003 to roug hly 200 MeV H in Eq. (|). 
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Also the bottom quark mass can be extracted from quantities related to the quark threshold. Recently pi|,pO[ 



precise value for the bottom quark mass has been determined in the context of QCD sum rules. For example, in |2C| ] 
the on-shell mass was eliminated in favour of the IS mass in order to reduce the error. Once M^^ is determined, 
mb{mi,) can be found in analogy with the top quark case. As a result the values M^^ = 4.71 ± 0.03 GeV and 
m\j{rrny) — 4.2 ± 0.06 have been obtained. Following the procedure described in j2^] one arrives at 

TOb(mb) = (4.71 - 0.40 - 0.11 - 0.03 ± 0.03 ± 0.04) GeV , (10) 

where the different terms correspond to different orders in the T-expansion. The first error is due to Ml^ and the 

second one reflects the error in a'p{Mz) ~ 0.118 ± 0.004 which is adopted from |2^. Due to the nice convergent 
behaviour in ( p^ the total error on mb(TOb) only contains these two sources which finally leads to 

TO6(m,,) = 4.17 ±0.05 GeV. (11) 

It is important to stress that taking into account of the newly computed 0{a'l) term in the MS-on-shell relation is 
crucial for the reliable estimation of the errors in (pi]). Indeed, a deviation of the real value for z'^ from the large- /3o 
estimation by, say, a factor of two, which one could not exclude a priori, would result to a systematic shift in mb(mf,) 
of around 100 MeV. 

A somewhat different approach for the determination of both the charm and bottom mass has been followed in . 
There the lower states in the heavy quarkonium spectrum were computed up to order and then used to extract the 
pole masses of h and c quarks. The transformation to the MS mass has been performed with the help of the two-loop 
relation However, to the order the quarkonium spectrum was computed it is more consistent to use the C'(a^) 
relation provided in this paper. Thus taking over the error estimates from [ p2[ the on-shell value for the bottom quark 
mass reads Mf, = ?>Mltollt GeV. It transforms to the following MS value 



mb(m6) = 4.322lO:^|3 Qgv, (12) 



where the value a^s\Mz) = 0.114 has been used as in |2^]. Compared to |22| the inclusion of the 0{a'l) terms leads 



to a shift in the central values of more than 100 MeV. In the case of the bottom quark this change is even larger than 
the errors presented in . This demonstrates that a consistent treatment of the different orders in Ug is absolutely 
crucial. 

To summarize: the computed value of the next-to-next-to-leading correction to the MS-on-shell mass relation proves 
to be in a good agreement with an estimation based on the large-/3o limit and has led to a significant reduction of the 
theoretical uncertainty in the determination of the quark masses. 

We would like to thank A.H. Hoang for useful conversation and I. Bigi, J.H. Kiihn and K. Melnikov for valuable 
suggestions and careful reading of the manuscript. This work was supported by DFG under Contract Ku 502/8- 
1 {DFG-Forschergruppe "Quantenfeldtheorie, C'omputeralgebra und Monte-Carlo-Simulationen") and the Schweizer 
Nationalfonds. 
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m(M)/M 




M/m(m) 






0{ai) 




o(«f) 


0{ai) 







-14.33 


-202(5) 


-11.67 


-170(5) 


13.44 


194(5) 


1 


-13.29 


-176(4) 


-10.62 


-146(4) 


12.40 


168(4) 


2 


-12.25 


-150(3) 


-9.58 


-123(3) 


11.36 


143(3) 


3 


-11.21 


-126(3) 


-8.54 


-101(3) 


10.32 


119(3) 


4 


-10.17 


-103(2) 


-7.50 


-81(2) 


9.28 


96(2) 


5 


-9.13 


-82(2) 


-6.46 


-62(2) 


8.24 


75(2) 



TABLE I. Dependence of Zm' and on n;. The choice fi^ = , respectively, fj? = has been adopted. ziS^ is defined 
as the sum of the terms inside the square brackets in Eq. (El). 



ni 


this work 


[23| (FAG) 


;2^ (PMS) 


[IC] (Iarge-/3o) 


2 


143(3) 


152.71 


153.76 


137.23 


3 


119(3) 


124.10 


124.89 


118.95 


4 


96(2) 


97.729 


98.259 


101.98 


5 


75(2) 


73.616 


73.903 


86.318 



TABLE II. Gomparison of the results obtained in this paper with estimates based of FAG, PMS and the large- /3o approxi- 
mation for M/m(m). 
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